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ABSTRACT 



Context. The PAMELA, Fermi and HESS experiments (PFH) have shown anomalous excesses in the cosmic positron 
and electron fluxes. A very exciting possibility is that those excesses are due to annihilating dark matter (DM). 
Aims. In this paper we calculate constraints on leptonically annihilating DM using observational data on diffuse ex- 
tragalactic 7-ray background and measurements of the optical depth to the last-scattering surface, and compare those 
with the PFH favored region in the moM — {(^av) plane. 

Methods. Having specified the detailed form of the energy input with PYTHIA Monte Carlo tools we solve the radiative 
transfer equation which allows us to determine the amount of energy being absorbed by the cosmic medium and also 
the amount left over for the diffuse gamma background. 

Results. We find that the constraints from the optical depth measurements are able to rule out the PFH favored region 
fully for the r~ -I- r"*" annihilation channel and almost fully for the /i~ -I- /i"*" annihilation channel. It turns out that 
those constraints are quite robust with almost no dependence on low redshift clustering boost. The constraints from 
the 7-ray background are sensitive to the assumed halo concentration model and, for the power law model, rule out 
the PFH favored region for all leptonic annihilation channels. We also find that it is possible to have models that fully 
ionize the Universe at low redshifts. However, those models produce too large free electron fractions at 2: > 100 and are 
in conflict with the optical depth measurements. Also, the magnitude of the annihilation cross-section in those cases is 
larger than suggested by the PFH data. 

Key words. Cosmology: theory - dark matter - diffuse radiation - cosmic microwave background - Elementary particles 



1. Introduction 

The PAMELA and Fermi satellite experiments and the 
HESS atmospheric Cherenkov telescope have shown an 
anomalous excesses in the cosmic electron and positron 
spectra. PAMELA has observed a steep rise of positron 
fraction e+/(e~-|-e+) at energies above 10 G eV with no sig- 
nifica nt excess in the cosmic antiproton flux (jAdriani et al.l . 
l2009t ). Fermi and HESS have measured an excess of 
high-energy (e~ + e"*") flux w ith a sharp cut-off of 
around 800 GcV (lAdriani et al.. 2009 yAbdo et all l2009t 
IHESS Collaboration: F. Aharonian . 200g). The ATIC and 
PPB -BETS balloon measurements i ndicat e a similar ex- 
cess (jChang et all . 120081: iTorii et al.l . l2008t) . The most ex- 
citing explanation to those anomalies is that they might be 
caused by the annihilation of dark matter (DM) particles. 
However, the nature of those signals require the proper- 
ties of DM to deviate strongly from the standard freeze-out 
predictions. The DM annihilation cross-section {<jav) has 
to be boosted some orders of magnitude over the standard 
freeze-out value of (cr^u)std — 3 x 10~^^cm'^/s and the an- 
nihilations should favorably occur only through the leptonic 
channels. 

Here we continue model independent studies of the 
DM annihilations into Standard Model charged leptons, 



DM + DM ^ £- + e+, e = e, /it, r, extending the ana l- 
jses pr esented in ser i es of works, e.g. ICirelli et al. (|2009[): 
BorricU o et all (I2009D: ICire Ui & Pali^ (|2009f ): iMeade et al] 
((2009i ): iArkani-Hamed et al.i (i2009|). Because DM annihila- 
tions to charged particles necessarily induce 7-ray signal, 
the observed 7-ray fluxes strongly constrain the DM anni- 
hilation scenario as an explanation to the above mentioned 
cosmic ray anomalies. Most of the 7-ray constraints arise 
from the observations of Galactic center where the DM den- 
sity is the highest. Those analyses take into account both 
the primary photons produced by final stat e radiation and 
decays of the DM an n ihilat ion products £ (jBertone et al.l 
120091 : iBeU k Jacauei I2009D as well as the secondary in- 
verse Compton (IC) scattering p hoton f lux produced by 
them (ICire lli fc Pa^ . 2009; Mead e et all [2009 ). 

In this study we focus on the possible effects of pho- 
tons produced in extragalactic DM annihilations all over 
the Universe today and in the past. We are particularly 
interested in finding out if the extragalactic photon sig- 
nals are able to rule out the PAMELA, Fermi, and HESS 
(P FH) favored r e gions on the m^M — {<^Av) plane as given 
in iMeade et al.l (|2009D . For this purpose we look at two 
extragalactic signals: (i) the diffuse 7-ray background in- 
cluding both the primary photons as well as the secondary 
photon flux produced in IC scattering of electrons off the 
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Cosmic Microwave Background (CMB) photons, (ii) the 
Thompson scattering optical depth of the CMB photons 
caused by the free electrons between us and the last scat- 
tering surface El- While the previous works on those top- 
ics take into account only the direct photons produced in 
DM annihilations, we show that in fact the diffuse IC pho- 
ton spectrum gives the strongest constraint. We stress the 
complementarity between those observables: the more en- 
ergy ends up in extragalactic gamma background the less 
energy is available for ionizing/heating the gas, and vice 
verse. The sp ectrum of diffuse 7-ray background has been 
measured by ISreekumar et al.l (119981 ) using data from the 
EGRET space experiment and the optical depth to decou- 
pling has been de rived from the CMB m easurements by the 
WMAP satellite (jPunklev et al.U2009[ ). New more precise 
data on 7-ray background is expected to be published by 
Fermi satellite this year. 

The effect of annihilating DM on reionization 
and recombination has b een previously investigated 
by a f ew authors, e. g. [ Padmanabhan fc Finkbeineil 



rriDM = 1 TeV 



I2OO5I); 'Mapelh et al.' ('2006'); 'Zhang et al." (200^; 
Nataraia n & Schwarz (2008); Bchkov & Hooper (2009_J; 
Galli et all (j20090 : ISlatver et al.l (I2009D. Prio r to th e 
calculations don e in iNataraian fc Schwarzl ([200I); 
iBelikov fc Hooperl (|2009a[ ) all the authors have only 
included the homogeneous DM component, i.e. have 
neglected the clustering effect. Despite the inclusion of the 
clustering there are some pr o blems with the formalism o f 
INataraian fc Schwarj ()2008D : IBelikov fc Hooped (|2009al) : 
the treatment for the e volution of the matter t emper ature 
is simplistic. Also, in INataraian fc Schward (|2008[ ). the 
optical depth f or the 7-rays is calculated incorrectly . The 
original work by iPadmanabhan fc Finkbeined (|2005|). and 
in par ticular, the extended newer version bv ISlatver et al.l 
(|2009r ) uses more accurate formalism, unfortunately the 
authors have not included the effect of DM clustering. 

During completion of our study couple of papers 
appeared calculating the diffuse 7-ray backgro u nd du e 
to the annihilating DM, e .g. iKawasaki et al.l (12009 ): 
Profum o fc Jeltemal (l2009l) : IBelikov fc Hooped (|2009b[ ). 
Kawas aki et al.l ~( 200E ) have neglected the 7-rays from 
the inverse Compton scattering of the final state e~ /e'^, 
which gives even str o nger bound than the primary 7-rays. 
iProfumo fc Jeltemal (|2009f) have used different formalism 
for the absorption of 7-ray s considering the numerical ap- 
proximation from IChen fc Kamionkowski (2004|). In this 
paper we include the relevant absorption processes directly 
(see the discussions below). 

Our paper is organized as follows: In Section 2 we de- 
scribe the source term, i.e. the energy input from the DM 
annihilation. In Section 3 we solve the radiative transfer 
equation and present a few example 7-ray spectra along 
with some results for the evolution of the ionization frac- 
tion and matter temperature. Section 4 presents our main 
results and our summary is given in Section 5. 



^ In the following the observational constraint derived from 
the Thompson scattering optical depth is called the "reioniza- 
tion constraint". 

^ The argument of exp in their Eq. (17). 





0.0012 - 




0.001 - 


LU 


0.0008 - 










T3 


0.0006 - 


LU 






0.0004 - 




0.0002 - 



electrons, DM + DM -- 
electrons, Df^ + DM -- 
electrons, DM + DM - 
photons, DM + DM 
photons, DM + DM 
photons, DM + DM - 
neutrinos, DM + DM - 
neutrinos, DM + DM - 



> e* + e", 96.0 

> + n", 34.0 
■ > t"*" + x", 1 5.9 
- > e^ + e", 4.0 
-- > + 2.8 




400 600 

E [GeV] 



1000 



Fig. 1. The final state distributions of electrons/positrons, 
photons, and neutrinos/antineutrinos for the annihilating 
DM with niDM = 1 TcV. 

Table 1. Energy partition between electrons/positrons, 
photons and neutrinos/antineutrinos for all three annihi- 
lation channels assuming rriDM = 100 GeV - 10 TeV. 



channel 


electrons 


photons 


neutrinos 


-1- e+ 


~ 96 -97% 


~ 3 - 4 % 


0% 




-34% 


~ 2 - 3 % 


~ 63 - 64 % 




~ 16% 


~ 16 - 17% 


~ 67 - 68 % 



2. Energy input from DM annihilation 

In our study we consider the DM annihilation to Standard 
Model charged lepton s following t h e lept onic annihilation 
scenario presented in ICirelli et al.l (|2009n . The final state 
distributions of electrons/positrons, photons, and neutri- 
nos/antineutrinos from the primary two-body states pro- 
duced in DM annihilations are co mput ed using MonteCarlo 
code PYTHIA (jSiostrand etTall . l2008t) . The example hnal 
particle distributions for all three annihilation channels as- 
suming the input DM particle with mass uidm = 1 TeV are 
shown in Fig. [T] The partition of energy between final par- 
ticles is shown in Table |T] The released electrons/positrons 
immediately interact with the CMB photons and upscatter 
those to high energies via the inverse Compton mechanism. 
Thus, in addition to the prompt photons released directly 
in the annihilation process, our photon spectrum contains 
the other part: the upscattered inverse Compton photons. 

Having an annihilating DM particle with mass itldm 
and with a thermally averaged cross-section {(Jav) the emit- 
ted power per volume and per steradian, i.e. the emission 
coefficient can be given as 



3{z 



1 {<yAv) 



(1) 



47r 2m£)M 

It is the emission coefficient for the homogeneous distribu- 
tion of DM. To include the effect of DM clustering p]jm{z) 

should be replaced by {p^Mi^)) = PDMiz)ii'^ + ^(^)) ) = 
B{z)p'j^]^_j{z). Here 6 = p/p — 1 is the density contrast 
and we have also defined a halo boost factor B{z) = 
{{l + 6{z))'^) = l + {S^{z)). To calculate the boost factor we 
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Fig. 2. (Upper panel) Boost factors B{z) as defined in 
Eq. ^) for various concentration models. (Lower panel) 
/-parameter as given by Eq. PH)) . i.e. the ratio of the to- 
tal energy deposition and local "smooth" energy injection 
rates for two different concentration models and two lep- 
tonic annihilation channels, assuming the annihilating DM 
particle with mass to_dm = ITeV. The short dashed lines 
show the ratio e/(47rj) where for calculating e one takes 
i? = 1 in Eq. ([5]). The results for the annihilation channel 
DM + DM ^ + T+ are very similar to the ^~ + 
case and for clarity arc not shown here. 



up at z ^ 100, which corresponds to the redshift where the 
DM halos with the assumed lowest masses start to form. 
As expected, the curves for the 10~^ Mq models start to 
deviate slightly earlier. It is clear that currently the biggest 
uncertainty influencing the magnitude of the boost factor 
is related to the way one chooses to extend the behavior of 
the concentration parameter below the mass scales directly 
probed by the simulations. We see that at low redshifts the 
variation between the models reaches two orders of magni- 
tude. On top of that there is an extra uncertainty related to 
the profiles of the DM halos. However, this is expected to 
be of smaller magnitude, and mainly due to those reasons, 
we have decided to use only NEW profiles throughout this 
paper. 

Using the boost factor we can write the total emission 
coefficient as 

j{z) = B{z)jiz). (4) 
We also define the corresponding spectral quantity 

Mz) ^ j{z)r{,y), (5) 

r(z/)dz. = 1, 

where the frequency distribution function f^{v) = + 
/prompt l'^) consists of two parts: (i) Inverse Compton part 
which is produced by the energetic electrons/positrons in- 
teracting wit h the CMB. Th i s is ca lculated following the 
formalism in ICirelli fc Pancil |2009'); (ii) Prompt photon 
part - these are the photons released directly from the an- 
nihilation event. 

For the inverse Compton scattering ffci^) = 
TTl/icr°(TTl) ^^'^ ^^"^ prompt emission /prompt (i^) = 



-0 
prompt 



use the halo model which approximates the matter distri- 
bution in the Universe as a superposition of DM halos (for 
a com prehensive review on halo model see iCoorav fc ShethI 
(|2002D ). Within this model B{z) can be given as 



B{z)^l 



A, 



3p„ 



oo 

- / dM^{M,z)f[c{M,z)] , 



(2) 



where Pm,o is the matter density at z = 0, = 200 is 
the overdensity at which the halos are defined, Mmin is the 
minimum halo mass, ■j^{M,z) is the halo mass function, 
c(M, z) represents the halo concentration parameter and 
the function /( c ) for t he halos with the NEW density profile 
(jNavarro et al.l . Il997t) IS given as 



m = 3 



1 



(1 



log(l + c) 



1 



(3) 



For the mass functio n we use the standard form given by 
ISheth fc TormenI (|1999f ). for the concentration parameter 
we use two models: (i) the model by Maccio et al. (2008), 
(ii) the power law model where c oc AI^^'^, which gives a 
good fit within the mass range resolved by the simulations. 
In our calculations we assume two values for the lower halo 
mass: 10~^ and 10~^ Mp, (fo r the motivation of those val- 
ues see e.g. lBringmannI ()2009D : Martinez ct al. (2009)). The 
resulting four models for the boost factor are given in the 
upper panel of Fig. O We see that the boost factor shoots 



3. Radiative transfer. Extragalactic gamma 
background. Reionization 

In the expanding Universe the formal solution of the ra- 
diative transfer equation for the intensity can be written in 
the form 



dz 

■ exp 



1 



1 + z 



(i + zY 



Hiz'){l + z') \l + z' 

Tu{z,Z ) 



dz^^±^B[z) 



(6) 



H{z') 



/ic(^)+/p™„pt(^) 



exp 



T^{Z,Z ) 



Here ii/{z) is the emission coefficient given by Eqs. ([I]), (U), 
(ini), H{z) is the Hubble parameter, factor oc (1 -I- z)~^ ac- 
counts for the cosmological dimming of the intensity, and 
the optical depth t^(2;, z ) is given as 



T^{z,z) = c i d 



H{z"){l + z")- 



(7) 



For similar comments on uncertainties related to the 

mass dependen ce of the concentration parameter see e.g. 

iMartinez et~aD l|2009f) . 
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()2005D . For their modification we need the total energy de- 
position rate per hydrogen particle which is given as 



£(z) = / eu{z)dv 




(9) 



To go beyond their "on the spot" approximation the /- 
parameter in their Eq. (5) s hould be replaced by the fol- 
lowing function of z (see also lSlatyer et alJ (|2009| i) 



I{z) 



<z) 
47rj(z) ' 



(10) 



1e-04 



0.01 



100 10000 
E^[mec2] 



1e+06 1e+08 



Fig. 3. The redshift z where the optical depth for photons 
reaches unity (i.e., Ti,{z,z ) = 1 in Eq.©) for several "ob- 
server's rcdshifts" : z = 0, 10, 100, 500, 1000. Here the energy 
plotted is the photon energy at redshift z. 



l+z 
1+z 



which gives us the ratio of the total energy deposition and 
local "smooth" energy injection rates. The function f{z) 
is plotted on the lower panel of Fig. [2] for two different 
concentration models, two leptonic annihilation channels 
{DM + DM er +e+, DM + DM fi^ +fj.+), assuming 
the annihilating DM particle with mass raoM = ITeV. The 
upturn of the curves at redshift z ~ 100 mimics the similar 
trend seen in the upper panel and is caused by the onset of 
the structure formation. The short dashed lines show the 
ratio 1/ (47rj) where for calculating e one takes i? = 1 in Eq. 
dni). These last two curves can b e directly coinpared with 
the results presented in Fig. 4 of ISlatver et all (|2009[ ). For 
clarity we have not shown the results for the annihilation 
channel DM + DM t~ +t'^ as those are very similar to 
v, and a^{z) is the absorption co- the -I- /Lt+ case. At large redshifts, where the Universe gets 

optically thick to photons, and where the structure boost 
B{z) — 1, one expects f{z) to asymptotically approach 
unity. This is indeed what we see in the case of + e+ 
annihilation channel. However, the asymptotic J{z) values 
for the /i^ + /i+ and t~ + r+ cases are smaller since a large 
fraction of energy is carried away by neutrinos (see Table 

In Fig. m we show some example 7-ray spectra at red- 
shift z = calculated with the above formalism for the 
annihilating DM particle with mass moM = ITeV. Here 
the two most extreme concentration models of Fig. [5] have 
been used and the results are given for all three leptonic 
annihilation channels considered in this paper. The ther- 
mally averaged annihilation cross-section has been set to 
25 times its standard value of ~ 3 x 10~^^cm^. The low 
energy bump of the characteristic double-bumped spec- 
trum is due to the inverse Compton process whereas the 
high energy bump is the prompt photon contribution. The 
points with errorbars correspond to the EGRET measure- 
ments of the ext r agalac tic gamma background as given in 
ISreekumar et al.l (|1998[ ). The solid horizontal line, which 
is used in the following as an upper bound for the level 
of diffuse 7-ray background, represents the EGRET mea- 
surements reduced by a factor of three to approximately 
account for the following: (i) dependent on the modeling 
details, it is predicted that from 25% up to 100% of the dif- 
fuse extragalacti c 7-ray bac k ground is due to u nresolved 
AGN s (.Stecker fc Salamonl . 1 19961 : IChiang fc M ukhcriei. 
119981: iMukheriee fc Chianei . 119991 : iMiicke fc Pohj i2000l) . 
(ii) the improved model for the G alactic contribution fur- 
ther suppresses the measurements ([Strong et al.ll200^ . (iii) 
there is an additional component due to the DM annihila- 
tion produced in our own Ga l actic halo (see the relevant 
calculation bv iKawasaki et all (|2009f) )FI 

* If the Reader has her/his own favored value for the possible 
reduction factor the final results in Fig. [6] are easily scalable to 



Here 

efficient. For av(z) we f ollow the formalism given in 
IZdziarski fc SvenssonI (jl989i) . The following processes are 
included in our calculations: 

— photoionization 

— Compton scattering 

— photon-matter pair production 

— photon-photon scattering 

— photon-photon pair production 

In Fig. [3] we plot the redshift z where the optical depth 
for photons reaches unity (i.e., t^(z,z ) = 1) for several 
"observer's redshifts": z = 0,10,100,500,1000. Here the 
energy plotted is the photon energy at redshift z. The low- 
est curve is the analog of t he curve plotted in Fig. 2 of 
IZdziarski fc Svensso i (119891) . 

Having a method to calculate the intensity of the photon 
field sourced by the annihilating DM at each redshift we can 
go on and express the energy deposition rate in the cosmic 
medium as 

e^{z) = AT:a^{z)h{z) . (8) 

To estimate the fraction of the deposited energy end- 
ing up in ionizing and heating the medium we use 
the approximation motiv a ted b y the original work 
of IShull fc van Steenberd (Il985f) and used in sev- 



eral subsequent papers (IChen fc KamionkowskiL 2004; 
Padman abhan fc Finkbcincr, 200 5|; iMapelh et all , boOC : 
Zhang et al.. . .2006 : .Nataraian fc Schwarzi 120081) : ^ (1 - 
Xe)/^ goes into ionization and (1 -I- 2xe)/'i into heating. 
Here Xe is the ionization fraction. To calculate the evolution 
of the ionization fraction and matter temperature we mod- 
ify the recombination code RECFAST (Scagcr ct al., 199' 
following the description in IPadmanabhan fc Finkbeine: 
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Fig. 4. Example 7-ray spectra at redshift z — assum- 
ing the annihilating DM particle with mass = ITeV. 
Here the two most extreme concentration models of Fig. 
[5] have been used and the results are given for all three 
leptonic annihilation channels. The thermally averaged an- 
nihilation cross-section has been set to 25 times its stan- 
dard value of ~ 3 X 10~^^cm^. The points with errorbars 
correspond to the EGRET measurements o f the extragalac - 
tic gamma background as given in SreckumareLall (| 19981) . 
The solid horizontal line, which is used in the following as 
an upper bound for the level of diffuse 7-ray background, 
represents the EGRET measurements reduced by a factor 
of three. 
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Fig. 5. The evolution of the ionization fraction (upper 
panel) and matter temperature (lower panel) for the model 
with the annihilating DM particle with mass moAi = ITcV 
and with a thermally averaged cross-section {aAv) = 3 x 
lO^^^cm^. The two most extreme concentration models of 
Fig. [2] have been used and for clarity only the results for 
the -|- channel are shown. The lowest solid line plots 
the case where the structure boost has been neglected, i.e. 
B[z) = 1 and the dashed line corresponds to the standard 
"concordance" cosmology without an annihilating DM. In 
the lower panel the dotted line shows the evolution of the 
CMB temperature. 



Note that due to relatively small fraction of energy re- 
leased in the form of prompt photons ^ 2 — 4%, except for 
the T~ + T"*" channel where the prompt photon part reaches 
up to 17%, the inverse Compton part of the spectrum has 
a higher amplitude than the prompt portion. As the obser- 
vational bounds are represented by a roughly horizontal line 
it is clear that in the case of e~ + and /i^ -I- ^.'^ channels 
the constraint we obtain on (cr^u) arises from the inverse 
Compton portion only. However, for the r~ -I- r"*" channel 
and for relatively low moMi when the prompt bump of the 
spectrum reaches the energy range probed by EGRET, the 
constraint on {uav) is provided by the amplitude of the 
prompt portion, instead. 

In Fig. [S] we show an example of the evolution of the 
ionization fraction (upper panel) and matter temperature 
(lower panel) for the model with the annihilating DM parti- 
cle with mass m^M = ITeV and with a thermally averaged 
cross-section {(Jav) = 3 x lO^^^'cm^. The two most extreme 
concentration models of Fig. [2] have been used and for clar- 
ity only the results for the + e+ channel are shown. 
The lowest solid line plots the case where the structure 
boost has been neglected, i.e. B{z) = 1 and the dashed 
line corresponds to the standard "concordance" cosmology 
without an annihilating DM. We have shown the evolu- 
tion only down to redshift z = 6 beyond which there is 
a clear observational evidence that the Universe gets fully 



accommodate that, since the DM annihilation signal is simply 
proportional to ((Taw). 



ionized (jFan et al.l . l2006l ). In the lower panel the dotted 
line shows the evolution of the CMB temperature. Note 
that at z ~ 100 the power law concentration model with 
Afmin = lO^^M© clearly starts to deviate from the smooth 
model without the clustering boost, resulting in an order 
of magnitude higher ionization fraction at low redshifts. 
However, if one calculates the optical depth to the last scat- 
tering then all of the models with DM annihilation plotted 
in Fig. [5] give very similar results, i.e. clustering boos t at 
low z has minimal effect. In lBelikov fc Hoopeij (|2009af) the 
authors discuss models where the DM annihilation is able 
to fully reionize the Universe at low redshifts. In our cal- 
culations we also find that, indeed, it is possible to fully 
reionize the Universe, however, all of those models give too 
large amplitudes for the ionization fraction plateau after 
z ~ 100 and thus violate the CMB measurements. 

4. Results 

In this section we are going to confront the leptonically an- 
nihilating DM models with the observational data. We are 
going to use: (i) measur ements of the extragala ctic 7-ray 
background by EGRET (|Sreekumar et al.l . ll998D . ( ii) mea- 
surement of the optical depth to the last scattering sur- 
face by WMAP (jDunklev et all . 120091 ). As described in the 
previous section ^ to account f or the better Galactic sub- 
traction due to IStrong et al.l (|2004[ ). for the annihilation 
signal from our own Galaxy, and for the contribution of 
AGNs to the 7-ray background, we use the upper bound of 



6 



Gert Hiitsi, Andi Hektor, Martti Raidal: Constraints on leptonically annihilating Dark Matter 



10000 




^<cluded by y background, power law c(M), M^■^^=^0 M^^l 
excluded by WMAPi 
excluded by Y background, Maccio et al. c(M), M^jn=10"^ 

PAMELA-Fermi-HESS favored region - 

GC-y- 

^ G-lcy- 



10000 - 



E 
o 



A 
> 

V 




1000 



100 



^excluded by ybackground, power law c(M), M^-^^= 

excluded by WMAP f 
excluded by ybackground, Maccio et al. c(M), M^^^=^0'^ 

PAMELA-Ferml-HESS favored region - 
GC-y- 

G-lcy- 



E 

o 



10000 



1000 



^ 100 



A 

> 



V 




100 



deluded by y background, power law c(r^), Mmjn="IO '^o^ 
excluded by WMAPl 
excluded by y background, Maccio et al. c(M), M|^jj,=10'^ M^ 
PAMELA-Fermi-HESS favored region - 
GC-y- 

G-lcy- 

' ' 1000 

rriDM [GeV] 



10000 



Fig. 6. Constraints from reionization and extragalactic 
gamma background for all three leptonic annihilation chan- 
nels. The upper, middle, and lower panel correspond to 
+ e+, /i^ + and + t+ channels, respectively. The 
solid line shows the PAMELA, Fermi, and HESS favored re- 
gion as given in Meade et al. (2009). We h ave also plotte d 
some Galactic constraints taken from Mea de et al.l (j2009t) : 
dashed lines show the bounds from the HESS measurements 
of the Galactic Center 7-rays, dotted lines are the bounds 
from FERMI measurements of 7-rays below 10 GeV. (See 
the main text for further details.) 



ISreekumar et al.l (|l998f ) reduced by a factor of three as our 
upper limit for the possible signal rising from the annihi- 
lating DM. This upper bound is shown as a solid horizontal 
hue in Fig. IH 



To calculate the constraint from reionization we use 
the WMAP measurement for the optical depth: r = 
0.087 ±0.017. Using the spectra of distant quasars we know 
that most of the hy drogen in the Uii iverse is ionized be- 
lo w redshift z ~ 6 (iFan et al.l . 20061). As was al s o done 
in lNataraian fc Schwar j (|2008D : iBelikov fc Hooper] (|2009i) 
we assume that helium is singly ionized below z — 6 and 
doubly ionized below z = 3. This assumption implies that 
the contribution to the WMAP optical depth from redshifts 
z > 6 is t{z > 6) = 0.047±0.017. In the following we use the 
1 — fT upper bound r(z > 6) = 0.064 for deriving our con- 
straints. The constraints we find on a moM — (ctav) plane 
for all three leptonic annihilation channels are given in Fig. 
[S] Here the upper, middle, and lower panel correspond to 
e~ +e^, /i" and +t+ channels, respectively. As it 

turns out, our constraints from WMAP r measurements are 
basically independent of the low redshift clustering boost, 
implying that already the high redshift (i.e. z > 100) en- 
ergy injection from the unclustered component is able to 
account for the measured level of t{z > 6). 

We also not e that, indeed, as pointed out in 
iBehkov fc Hooper! ^QQM) one can have models where the 
annihilating DM is able to fully reionize the Universe at low 
redshifts. However, we find that those models lead to the 
strong violation of the WMAP optical depth measurement 
due to the vastly increased level of the ionization fraction 
plateau beyond redshifts z ^ 100. 

In Fig. [6] we have also plotted region s favored by the 
PFH data as given in i Meade et all (|2009l ). We see that for 
the power law concentration model with -Mmin = 10~^ Mq 
and for all three annihilation channels the constraints from 
the extragalactic 7-ray background completely exclude the 
PFH favored region. This is also the case for the power law 
c{M) model with Mmin = 10~^ Mq, which for clarity is 
not shown in the figure. The WMAP r constraint, which is 
insensitive to the structure boost, and thus does not depend 
on large uncertainties rela ted to c(-M), i s able to rule out 
the PFH favored region of iMeade et al.l (|2009l ) in the case 
of T~ + r+ channel and also by a large part in the case 
of M~ + M"*" a nnihilation channel. The excluded part of the 
iMeade et al.l (2009.) favored region in the e~ -I- e"*" case is 
somewhat smaller. 

In addition to the extragalactic constraints in Fig. [6] we 
have also shown tw o types of Galactic bounds as found in 
IMeade et all (|2009D : (i) constraints derived from the HESS 
observations of the Galactic Center 7-rays (0.1° angular 
region), plotted by dashed lines and labeled as GC-7, (ii) 
constraints from the FERMI observations of the 7-rays be- 
low 10 GeV in the 10° < |6| < 20° region (i.e., away from 
the Galactic Center), plotted by dotted lines and labeled 
as G-IC70 Although both of those Galactic bounds shown 
here assumed NFW density profile, the G-IC7 bound is 
less dependent on this assumption, as the observed region 
is away from the Galactic Center. It is evident that in sev- 
eral occasions the complementary extragalactic constraints 
derived in this work provide stronger bounds on annihi- 
lation cross-section, and thus are certainly of considerable 
interest. 

As a final point, note that our constraints based on dif- 
fuse extragalactic 7-ray background depend significantly on 
the assumed fraction of unresolved astrophysical sources 



At those lower energies the Inverse Compton part of the 
annihilation 7-spectrum dominates. 
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contributing to the EGRET signal. There is a general 
consensus that a significant fraction of the diffuse signal 
is due to unresolved AGNs, however the predicted val- 
ues are highly model dependent with the fractions rang- 
ing from 25% up to 100%. The upcoming measurements 
by Fermi satellite should clarify this situation considerably, 
since compared to EGRET, Fermi has significantly bet- 
ter sensitivity and angular resolution, allowi ng it to resolve 
many more individual astrophysical sources (jAtwood et al.l . 
120091) ■ 

5. Summary 

In this paper we have constrained observati onally moti- 
vated leptonic DM annihilation scenario of ICirelli et al.l 
(|2009D using measurements of the extrag alactic 7-ray back- 
groun d by EGRET space experiment (jSreekumar et al.l . 
I1998 D and the measurement of the optical depth to decou- 
phng by WMAP satellite (Dunklcy et al., 2009). Our main 
results for all considered leptonic annihilation channels are 
given on three panels of Fig. [5] where we ha ve also added 
the re gion favored by the PFH data as given in lMeade et al.l 

iooi). 

Here are our main conclusions: 

— For the power law concentration model c oc M^'^'^ the 
constraint from the extragalactic gamma background 
rules out the PFH favored region of lMeade et all (|2009[ ) 
for all three annihilation channels considered in this 
work. 

— The same constraint in the case of[M accio et al.l (I2008D 
concentration model is significantly weaker owing to the 
much smaller structure formation boost as seen on the 
upper panel of Fig. [51 Compared to the mass of the 
smallest DM halos assumed in this work numerical sim- 
ulations probe several orders of magnitude larger mass 
scales. The particular way one chooses to extrapolate 
the c(M) relation to smaller scales is currently arguably 
the biggest source of uncertainty for the annihilation- 
sourced 7-ray background calculations. 

— For the minimal DM halo masses assumed in this pa- 
per the reionization constraints, on the other hand, turn 
out to be insensitive to the low redshift clustering boost, 
and thus are free of this large uncertainty in the mass- 
concentration relation. Thus these constraints are sig- 
nificantly more robust. 

— The primary photons give important contribution to the 
derived constraints only in the case of + r+ annihi- 
lation channel. In the case of /i~ -I- /i+ and + e+ 
annihilation channels the constraints are mostly given 
by the IC scattering photon flux. 

— For the + annihilation channel the constraints 
from reionization completely rule out the PFH favored 
region whereas for the fi~ + case large part of the 
region gets excluded. In the case of e~ + annihilation 
channel the excluded part is relatively small, but one 
has to keep in mind that in reality this channel does 
not provide a part icularly good fit to the PFH data 
(|Meade et al.l . l2009D. 

— In agreement to lBelikov fc Hoopeil (l2009aD we find that, 
indeed, it is possible to have annihilating DM models 
that completely ionize the Universe at low redshifts. 
However, we find that those models result in a too high 
ionization fraction plateau at z > 100, leading to the 



optical depths which are incompatible with the WMAP 
measurements. Also, the annihilation cross-section one 
needs is somewhat larger than suggested by the PFH 
data. 

— In the light of our results, the foreseen new diffuse 7-ray 
data by Fermi satellite is expected to further constrain 
the DM annihilation solutions to the cosmic ray anoma- 
Hes. 
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